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A C# package is presented that allows a user for an input quantum circuit to generate a set of multivariate 
polynomials over the finite field Z2 whose total number of solutions in Z2 determines the output of the quantum 
computation defined by the circuit. The generated polynomial system can further be converted to the canonical 
Grobner basis form which provides a universal algorithmic tool for counting the number of common roots of the 
polynomials. 



1. INTRODUCTION 

One important aspect of quantum computation 
is estimation of computational power of quantum 
logical circuits. As it was recently shown in pQ, 
determining the output of a quantum computa- 
tion is equivalent to counting the number of solu- 
tions of a certain set of polynomials defined over 
the finite field Z2. 

Using ideas published in £Q, we have written 
a C| program enabling one to assemble an ar- 
bitrary quantum circuit in a particular universal 
gate basis and to construct the corresponding set 
of polynomial equations over Z2 . The number of 
solutions of the set defines the matrix elements of 
the circuit and therefore its output value for any 
input value. 

The generated polynomial system can further 
be converted into the canonical Grobner basis 
form by applying efficient involutivc algorithms 
described in 0. A triangular Grobner basis for 
the pure lexicographical order on the polynomial 
variables is generally most appropriate for count- 
ing the number of common roots of the polyno- 
mials. 

Our program has a user-friendly graphical in- 
terface and a built-in base of the elementary gates 
representing certain quantum gates and wires. A 
user can easily assemble an input circuit from 
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those elements. 

The structure of the paper is as follows. In 
Section [5] we outline shortly the circuit model of 
quantum computation. Section presents the fa- 
mous Feynman's sum-over-paths method applied 
to quantum circuits. In Section 0] we describe a 
circuit decomposition in terms of the elementary 
gates. In Section [5] we show how to assemble an 
arbitrary circuit composed from the Hadamard 
and Toffoli gates that form a universal basis. Sec- 
tion[S]demonstrates a simple example of handling 
the polynomials associated with a quantum cir- 
cuit by constructing their Grobner basis. We con- 
clude in Section [7| 

2. QUANTUM CIRCUITS 

To quantize the classical bit, we go from the 
two-element set {0, 1} to a two-level quantum 
system described by the two-dimensional Hilbert 
space C 2 . In contrast to the classical case, the 
quantum bit (qubit) can be found in a superpo- 
sition \ip) = a\0) + f3\l) of the states |0) and |1) 
called a computational basis, where a, (3 S C are 
the probability amplitudes of |0) and |1) respec- 
tively. 

The simplest quantum computation is a uni- 
tary transformation on the qubit state 

\<p) = U\il>), Urf = I. 

A measurement of the qubit in the computational 
basis |0) and |1) transforms its state to one of the 
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basis states with probabilities determined by the 
amplitudes 



a 



|0)+/3|1) 



|0) with probability \a\ 2 
|1) with probability |/3| 2 



To compute 
function / : Z 2 l 



a reversible Boolean vector- 
-> Z2, one applies the appropri- 
ate unitary transformation Uf to an input state 
I a) composed of some number of qubits 



Uf |a) , 



lb) G 



The output state |b) is not the outcome of the 
computation until its measurement. After that 
the output state can be used anywhere. 

Some unitary transformations are called quan- 
tum gates. A quantum gate acts only on a few 
qubits, on the rest it acts as the identity. A quan- 
tum circuit can be assembled by appropriately 
aligning quantum gates. The unitary transfor- 
mation defined by the circuit is the composition 
of the constituent unitary transformations 



U f = U m U m - 1 ---U 2 U 1 



(1) 



A quantum gate basis is a set of universal quan- 
tum gates, i.e. any unitary transformation can 
be presented as a composition of the gates of the 
basis. As well as in the classical case, there are 
several sets of universal quantum gates. For our 
work it is convenient to choose the particular uni- 
versal gate basis consisting of Hadamard and Tof- 
foli gates 

The Hadamard gate is a one-qubit gate. It 
turns a computational basis state into the equally 
weighted superposition 
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The resulting superpositions for |0) and |1) differ 
by a phase factor. 

The Toffoli gate is a tree-qubit gate. Input bits 
x and y control the behavior of bit z, and the 
Toffoli gate acts on computational basis states as 

(x,y,z) i-> (x,y,z®xy) 

An action of a quantum circuit can be de- 
scribed by a square unitary matrix whose matrix 



element (b|[//|a) yields the probability ampli- 
tude for transition from an initial quantum state 
I a) to the final quantum state |b). The matrix 
element is decomposed in accordance to the gate 
decomposition of the circuit unitary transforma- 
tion <£Q) and can be calculated as sum over all the 
intermediate states a;, i = 1,2, . . . m - 1: 
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b| U m |a„ 



(ailEMa) 



3. FEYNMAN'S SUM-OVER-PATHS 

To apply the famous Feynman's sum-over- 
paths approach to calculate the matrix element 
of a quantum circuit, we replace every quantum 
gate of the circuit under consideration by its clas- 
sical counterpart. The trick here is to select 
the corresponding classical gate for the quantum 
Hadamard gate because for any input value, or 
1, it gives with equal probability either or 1. We 
denote the output of the classical Hadamard gate 
by the path variable x. Its value determines one 
of the two possible paths of computation. The 
classical Toffoli gate acts as 

(di, a 2 , a 3 ) i-> (01, a 2 , a 3 © a x a 2 ) , 

and the classical Hadamard gate as 

a\ 1— > x a,i, x G Z2 



— H 

— H 




H 










-€ 


3— 


r 


V- 


H 













a, 
a 1 



— H 

— H 




H 




.r, © x,x 4 








X, 




a, 8 x,x 2 


H 













bJx) 



U- b 3 (x) 



x = (x, . x, . x 3 , x 4 ) r e Zi 



Figure 1. From quantum to classical circuit 
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Fig. n shows an example of quantum circuit 
(taken from PP) and its classical correspondence. 
The path variables Xi comprise the (vector) path 
x = (xi,X2, X3, x A ) T G Z|. 

A classical path is a sequence of classical bit 
strings a, a 1, 02, ... , a rn = b resulting from appli- 
cation of the classical gates. For each selection 
of values for the path variables xi we have a se- 
quence of classical bit strings which is called an 
admissible classical path. Each admissible classi- 
cal path has a phase which is determined by the 
Hadamard gates applied. The phase is changed 
only when the input and output of the Hadamard 
gate are simultaneously equal to 1, and this gives 
the folmula 

ip{x) = input • output 

Hadamardgatcs 

Toffoli gates do not change the phase. 

For our example the phase of the path x is 

<ys(x) = aixi © a 2 X2 ffi X1X3 © 2:4(03 © 21X2) 

The matrix element of a quantum circuit is 
given by sum over all the allowed paths from the 
classical states a to b 



(b| U f |a) = 



1 



E (-D 



ip(x) 



x:b(x)=b 



where h is the number of Hadamard gates. The 
terms in the sum have the same absolute value 
but vary in sign. 

Let No be the number of positive terms in the 
sum and Ni the number of negative terms 
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\{x\b(x) 
\{x\b(x) 



& (p(x) 
& (p{x) 
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These equations count solutions to a system of 
n+1 polynomials in h variables over Then the 
matrix element may be written as the difference 



(b| U f |a) = 



1 



(No - JVi) 



4. CIRCUIT DECOMPOSITION 

To provide a user with a tool for assembling 
arbitrary quantum circuits composed from the 
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Figure 2. Circuit decomposition into elementary 
gates 



Hadamard and Toffoli gates we represent a cir- 
cuit as a rectangular table (Fig.|2J|. 

Each cell in the table contains an elementary 
gate from following set 



E 



+ A v A V A V 

{I,I,I,I,M,M,A,A,H} 



(2) 



so that the output for each row is determined by 
the composition of the elementary gates in the 
row. Thereby, each elementary unitary transfor- 
mation Uj is represented as an n-tuple of elemen- 
tary gates. 

Fig. shows action of the elementary gates 
from J2J): the identities, the multiplications, the 
additions modulo 2, and the classical Hadamard 
gate. The identity just reproduces its input. The 
identity-cross reproduces also its vertical input 
from the top elementary gate to the bottom one 
and vice versa. Every identity-down and identity- 
up have two outputs - horizontal and vertical. 
The multiplication-up and multiplication-down 
perform multiplication of their horizontal and 
the corresponding vertical inputs. In a similar 
manner act the addition-up and addition-down. 
Each Hadamard gate outputs an independent 
path variable irrespective of its input and can give 
a nonzero contribution to the phase. 
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phase is also calculated. Then we proceed the 
same way with the second column, with the third 
column, and so on up to the last column. Fig. 
shows an example. 
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Figure 4. Assembling circuit 
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Figure 3. Action of elementary gates 



5. ASSEMBLING CIRCUITS 

How can one assemble a circuit? First of all, 
we define an empty table of the required size. In 
this case both output and phase are not fixed. 
Then we place the required elementary gates in 
appropriate cells. Now the output is the result of 
applying the elementary gates to the input. The 



Apart from the ordinary menu, our program 
contains the toolbar for selecting an elementary 
gate and the toolbar for main operations. There 
are two windows: for assembling a circuit and for 
showing its polynomials. 

A circuit is represented in the program as 
two 2c?-arrays: one for the elementary gates and 
another for their polynomials. The phase poly- 
nomial is separately represented. The following 
piece of code demonstrates construction of the 
circuit polynomials 

for each Column in Table of Gates 
for each Gate in Column { 
construct Gate Polynomial ; 
if Gate id Hadamard 

reconstruct Phase Polynomial; } 

The method for constructing a gate polynomial 
is recursive because of the need to go up or down 
for some gates. 

Any circuit is saved as two files. One file is bi- 
nary and contains the circuit itself. Another file 
has a text format. It contains the circuit poly- 
nomials in a symbolic form. The program allows 
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to save polynomials in several formats convenient 
for loading into a computer algebra system (for 
example, in Maple or Mathematica) for the fur- 
ther processing. It is also possible to load back in 
memory a saved circuit. 

Note, that the part of our code for name space 
PolynomialJVIodulo_2 written in C# 4 can 
also be used independently on our program. This 
part contains classes for handling polynomials 
over the finite field Z2. Class Polynomial is a list 
of monomials, class Monomial is a list of letters, 
class Letter is an indexed letter provided with a 
positive integer superscript (power degree). 

6. QUANTUM POLYNOMIALS 

A system generated by the program is a finite 
set F C R of polynomials in the ring 

R := Z 2 [a. l ,b j ][xi, ...,x h ] 
aj, bj e Z 2 , i, j — 1, ...n 

in h variables and 2n binary coefficients. One has 
to count the number of roots No and Ni in Z2 of 
the polynomial sets 

F = {f,...,f k ,<p}, F 1 ={/,... + 1} 

Then the circuit matrix is given by 

1 



(b| U I a) = 



(JV0-JV1) 



To count the number of roots one can convert 
Fq and F\ into a triangular form by computing 
the lexicographical Grobner basis by means of the 
Buchberger algorithm or by involutive algorithm 
decribed in 

For the example shown on Fig. ^ we have the 
following polynomial system: 

fi = x 2 Xi + x 3 + bi 
f2 = x 2 +b 2 
f3 = X4 + b 3 

(p = x x x 2 + xix 3 + a%x\ + a 2 x 2 + a 3 x 4 

The lexicographical Grobner basis for the order- 
ing x\ >- x 2 y x 3 >- x 4 on the variables and 
representing both Fo and F\ is as follows 

91 = {a\ + h)xi + a 2 b 2 + a 3 b 3 (+1) 

92 =x 2 +b 2 

.93 = x 3 + bi+ b 2 b 3 
g 3 = x 4 + b 3 



From this lexicographical Grobner basis we im- 
mediately obtain the following conditions on the 
parameters: 

<2i + 61 = & a 2 b 2 + a 3 b 3 = 
a\ + bi = & a 2 &2 + 03^3 = 1 

From these conditions we easily count 2 (0) roots 
of F (Fx) and (2) roots of F (Fx). In all other 
cases there is 1 root of Fo and Fi . 
Some matrix elements are 

(000| U |001) = ^ , (000|Z7|111) = 

7. CONCLUSION 

We presented the first version of a program tool 
for assembling arbitrary quantum circuits and for 
constructing the corresponding polynomial equa- 
tion systems. Its number of solutions uniquely 
determines the circuit matrix. 

There is the algorithmic Grobner basis ap- 
proach to converting the system of quantum poly- 
nomials into a triangular form which is useful for 
computing the number of solutions. 

Thus, the above presented software together 
with Grobner bases provide a tool for simulating 
quantum circuits. 
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